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OPTICAL  METHOD  FOR  MEASUREMENT  OF  ULTRASONIC  REFLEC¬ 
TION  FROM  SOLID  PLATES. 

Walter  G.  Mayer  and  Tran  D.K.  Ngoc. 

Department  of  Physics,  Georgetown  University,  Washington,  DC  20057,  USA. 

Abstract 

An  acousto-optic  interaction  method  was  used  to  measure  the  modulus  of  the  reflec¬ 
tion  for  ultrasonic  waves  impinging  at  various  angles  on  solid  plates  immersed  in 
weter.  Theor\  and  experiment  are  shown  to  agree  qualitatively  when  the  product 
frequency  times  plate  thickness  is  small. 

Introduction 

The  reflection  coefficient  of  ultrasonic  waves  from  a  solid  plate  immersed  in  a  liquid 
has  been  given  by  Pitts  et  al.  ( 1  ]  for  the  simplified  case  where  absorption  is  neglected. 
The  general  form  of  this  amplitude  reflection  coefficient 

RA=N/FaFs  (1) 

is  the  same  when  absorption  is  included,  as  was  done  by  Ngoc  12] ,  in  which  case  the 
quantities  in  eq. (1)  are  defined  in  the  following  manner. 


N  =  N?  +  N| 


+  2N|Nj(l  -  cosPcosQ)/(sinPsinQ), 


Fa  =  ((1  -cosP),'sinP]  N|  +  [( 1  -cosQ)/sinQ]  N2  +iNj,  (3) 

Fs  =  |(1  +  co$P),sinP]  N,  +  1(1  +  cosQl.'sinQ]  N2  -  iNa,  (4) 

N,=  (kl  -  2k-)2,  (5) 

N2=4kiKdKs.  (6) 

N3=pk*KdK,  (7) 

with  the  sy  mbols  in  eqs.  (2)-(7j  gj\en  by 
k  =  (  Ztti  \  1  ( 1  +  ia/2rr).  \ 

kd  -  f  - rf  vd)  (1  +  iad  Z~).  I  (g) 

ks  =  Crf  vs)  (i  +  jas  Zr).  ( 

kx  =  k  sing,.  I 

K  =  Ik"  k^ )  ", 

Kd  =  (ka  -  k-)‘:.  (9) 

Ks  -  (k^  kx I  , 

P  =  dKd  I  (10) 

0  =  dks  j 

Here  d  is  the  thickness  of  the  solid  plate  (in  mm)  and  f  the  ultrasonic  frequency 
(in  MHr).  The  angle  of  incidence  is  denoted  by  0,,  and  p  is  the  ratio  of  the  liquid 
density  to  the  solid  density.  The  parameters  v.  vd,  vs.  and  a,  ad.  as  are  the  sound 
velocities  and  attenuation  per  wavelength  of  the  wave  in  the  liquid,  and  the  longitud¬ 
inal  and  shear  waves  in  the  solid,  respectively. 

Numerical  evaluation  of  the  reflection  coefficient,  eq.  ( 1 ).  shows  that  reflec¬ 
tion  peaks  occur  for  some  0,  with  the  number  of  peaks  generally  increasing  when  the 
product  fd  increases.  The  height  of  these  peaks  and  the  minima  of  RA  between  peaks 
depend  strongly  on  the  values  of  ad  and  a^.  The  reflectivity  pattern  is  rather  compli¬ 
cated  when  fd  becomes  larger  than  about  10,  making  the  experimental  measurement 
of  Ra  as  a  function  of  8t  difficult  due  to  the  existence  of  nor.-specular  reflection  1 3 1 . 

However,  when  thin  plates  are  selected  ihe  results  of  calculations  based  on 
eq  (  I )  can  be  ver.fied  experimentally  using  an  acousto-optic  interaction  technique. 
This  method  and  the  results  are  described  below  . 
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Figure  I  (a)  Schematic  diagram  of  beam  alignment  w  ith  respect  to  solid  reflector  (S). 

Incident,  reflected,  and  transmitted  beams  are  denoted  by  J,  R,  and  T,  res¬ 
pectively. 

(b)  Schematic  diagram  of  alignment  of  first  and  second  orders  of  diffrac¬ 
tion  patterns  produced  by  incident  and  reflected  ultrasonic  beams  shown 
in  (a) 


Method  and  Results 

The  experimental  setup  consists  of  a  basic  Raman-Nath  [4]  diffraction  arrangement, 
part  of  which  is  shown  schematically  in  Fig.  la.  A  laser  beam  is  expanded  into  a 
collimated  beam  of  approximately  10  cm  diameter,  indicated  by  the  circle  in  Fig.  la. 
This  beam  travels  through  a  glass  tank  filled  with  water,  and  the  upper  surface  of  the 
solid  plate  reflector  tS)  is  located  in  the  center  of  the  beam.  This  surface  is  parallel 
to  the  propagation  vector  of  the  light  so  that  an  ultrasonic  incident  beam  (I),  the 
reflected  beam  (R)  and  the  transmitted  beam  (T)  are  at  right  angles  to  the  light  prop¬ 
agation  direction,  satisfying  the  simplest  conditions  of  light  diffraction  by  ultrasonic 
waves  as  given  by  Raman  and  Nath  [4] ; 

A  large  lens  (about  10  cm  diameter)  is  used  to  focus  the  light  beam  after  it 
has  traversed  the  water  in  the  glass  tank.  With  no  sound  field  present  in  the  path 
of  all  the  light  will  be  focused  in  one  point,  indicated  by  0  in  Fig.  lb.  How¬ 
ever,  if  1,  R.  and  T  are  present,  diffraction  patterns  will  be  formed  centered  around  0. 
The  alignment  of  the  diffraction  pattern  caused  by  the  incident  beam,  indicated  by 
21.  II,  0,  etc.  depends  on  the  angle  of  incidence  of  the  beam  f,  and  the  alignment  of 
the  pattern  created  by  the  reflected  beam  (2R.  1  R,  0,  etc)  on  the  direction  of  travel 
of  R.  The  transmitted  beam  T  will  produce  a  diffraction  pattern  w  hich  is  super¬ 
imposed  on  the  pattern  generated  by  1,  since  1  and  T  have  the  same  direction  of 
travel.  Thus,  if  one  blocks  all  the  light  traversing  the  water  below  the  solid  plate, 
indi.ated  by  the  dashed  line  in  Fig  1  a,  one  obtains  a  set  of  diffraction  r3tterns 
caused  only  by  I  and  R. 

The  light  intensity  contained  in  the  first  orders  of  the  diffraction  pattern, 
i  e. ,  I  ]  |  or  I  j  r,  depends  on  the  amplitudes  of  the  two  ultrasonic  beams  I  and  R. 

Since  the  sound  amplitude  of  the  reflected  beam  changes  as  a  function  of  the  angle 
of  incidence,  as  predicted  from  eq.  ( 1 ).  the  Light  intensity  1  j  r  also  changes  accord¬ 
ing  to  the  theory  of  Raman  and  Nath.  This  light  intensity  is  given  by 

Jir  =  J?('r),  01) 

where  Jt  is  the  first-order  Bessel  function,  whose  argument  vr  is  proportional  to 
the  amplitude  of  the  reflected  ultrasonic  beam. 

If  one  now  keeps  vj,  the  amplitude  of  the  incident  beam,  constant  as  one 
changes  the  angle  of  incidence,  one  can  compare  the  changing  Ijr  to  the  fixed  I  j  j . 
One  can  thus  plot  the  amplitude  of  the  reflected  ultrasonic  beam  in  comparison 
to  the  incident  amplitude.  This  procedure  yields  a  plot  of  reflectivity  for  a  given 
solid  reflector  as  a  function  of  6,. 
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Angle  of  Incidence 


Figure  2  Calculated  values  (dashed  line)  and  experimentally  measured  values 
(circles)  of  reflection  coefficient  for  a  Plexiglas  plate  in  water.  Freq¬ 
uency  2  MHz.  plate  thickness  1/16  inch,  fd  =  3.1  75  MHz. mm 


Angle  of  Incidence 


Figure  3  Calculated  values  (dashed  line)  and  experimental!)  measured  values 
(circles)  of  reflection  coefficient  for  a  Plexiglas  plate  in  water.  Freq¬ 
uency  2  MHz,  plate  thickness  1  mm,  fd  =  2  MHz. mm 

The  results  for  two  thin  Plexiglass  plates  in  water  are  shewn  in  Figs  2  and  3. 
The  dashed  lines  are  the  theoretically  calculated  curves  for  the  reflection  coefficient 
where  the  computation  is  based  on  eq.  (1).  Qualitative  agreement  exists,  partic¬ 
ularly  referring  to  location  and  number  of  reflection  maxima  and  minima  The 
quantitative  agreement  is  marginal,  caused  primarily  by  uncertainties  about  the 
exact  value  of  the  absorption  coefficient  for  both  the  longitudinal  and  shear  wave 
absorption  coefficient  of  Plexiglas.  An  analysis  of  eq  (1 )  shows  that  the  heights 
of  the  reflection  coefficient  peaks  and  the  depth  of  the  minima  are  extremely 
strong  functions  of  the  absorptivity,  while  the  location  and  the  number  of  extrema 
are  not  materially  influenced  by  the  value  of  the  absorptivity. 

It  should  be  noted  that  the  method  described  here  is  applicable  when  the 
reflected  beam  has  an  amplitude  cross-section  which  is  essentially  the  same  as  that 
of  the  incident  beam.  This  condition  is  satisfied  in  most  cases  where  fd  is  small 
However,  when  fd  becomes  large,  non-specular  reflections  and  phase  shifts  between 
portions  of  the  reflected  beam  may  occur  (5).  and  the  simple  form  of  the  uniform- 
beam  Raman-Nath  theory  may  no  longer  hold 
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Conclusion 

The  acousto-optic  interaction  method  described  above  can  be  used  to  determine  the 
locations  of  reflection  coefficient  maxima  and  minima  for  ultrasonic  beams  imping¬ 
ing  at  various  angles  on  solid  plates.  The  method  yields  qualitative  results  of  the 
changes  in  the  amplitude  reflection  coefficient  of  plates  in  a  liquid  provided  phase 
shifts  within  the  reflected  beam  are  not  severe. 

Acknowledgement 

This  work  was  supported  by  the  Office  of  Naval  Research,  U.S.  Navy. 

References 

11]  L.E.  Pitts,  T.J.  Plona.W  G.  Mayer.  lEEE-Trans.  Sonics  &  Ultrs.  SU-24,  101 
(1977).  -  ’ 

[2]  T  D  K.  Ngoc,  Ph.D.  Thesis.  Physics  Department,  Georgetown  University 
Washington,  DC,  (1979). 

13]  M.  Behravesh,  L.E.  Pitts,  W.G.  Mayer.  Acoustics  Letters,  1,  165,(1977). 

[4]  C.V.  Raman  and  N.S.  Nath.  Proc.  Indian  Acad  Sci.  A2,  406,(1935). 

[5]  W.G.  Neubauer.  J.  Appl.  Phys.  44,  48,(1973). 

(Received  4  June  1980) 


Aiousncs  Leaeri  I'ol  3.  \o  9,  1980 


IEEE  TRANSACTIONS  ON  SONICS  AND  ULTRASONICS,  VOL.  SU-27,  NO.  S,  SEPTEMBER  1V80 


229 


A  General  Description  of  Ultrasonic  Nonspecular 
Reflection  and  Transmission  Effects  for 

Layered  Media 

TRAN  D.  K.  NGOC  and  WALTER  G.  MAYER 


Abstract- A  numerical  integration  method  is  used  to  calculate  the  in¬ 
tensity  distribution  in  reflected  and  transmitted  beams  for  liquid-solid- 
liquid  layered  media.  It  is  shown  that  the  formulation  describes  all 
nonspecular  reflection  and  transmission  phenomena  for  all  angles  of 
incidence,  including  critical  angles.  It  is  found  that  the  existence  of 
nonspecular  phenomena  not  only  depend  on  the  angle  of  incidence  and 
the  product  frequency  times  solid  thickness  but  also  on  the  product 
frequency  times  beamwidth.  Various  examples  are  given  for  reflection 
and  transmission  for  incident  angles  corresponding  to  critical  angles, 
plate-mode  angles,  and  be  tween-mode  angles. 

I.  Introduction 

EOMETRICAL  reflection  or  transmission  of  an  ultra¬ 
sonic  bounded  beam  produces  a  beam  with  probably  a 
lower  intensity  level  than  the  incident  beam  but  with  the  same 
intensity  distribution  profile.  This  is  called  specular  reflectivity 
or  transmittivity.  Nonspecular  reflection  or  transmission  ef¬ 
fects  refer  to  phenomena  where  the  reflected  or  transmitted 
beam  has  an  intensity  distribution  profile  different  than  that 
of  the  incident  beam;  this  may  include  a  lateral  beam  displace¬ 
ment,  one  or  several  minimum  intensity  areas,  and  a  trailing 
sound  field.  Fig.  1  illustrates  some  typical  profiles  of  a  beam 
nonspecularly  reflected  from  or  transmitted  through  a  layered 
system. 

Bertoni  and  Tamir  [1  ]  successfully  developed  an  analytical 
model  to  describe  nonspecular  reflectivity  for  a  liquid-solid 
system  at  Rayleigh  angle  incidence.  Later  Bertoni  and  Hou 
(2]  slightly  modified  this  analysis  to  incorporate  attenuation 
in  the  media.  In  the  same  analytical  framework,  Breazeale  et 
al.  (3)  examined  the  reflected  profile  observed  at  large  dis¬ 
tances  from  the  liquid-solid  interface.  Pitts  et  al.  [4]  ex¬ 
tended  Bertoni  and  Tamir’s  model  to  investigate  the  profile  of 
a  bounded  beam  reflected  from  a  solid  plate  immersed  in  a 
liquid.  The  same  approach  was  also  used  by  Ng  [5]  to  study 
the  transmitted  profile  for  the  same  layered  system.  Their  cal¬ 
culated  results  are  qualitatively  consistent  with  experimental 
observations  but  only  applicable  when  the  resonant  conditions 
prevail,  that  is,  for  Rayleigh  angle  incidence  in  the  liquid-solid 
case  and  for  plate-mode  incidence  in  the  liquid-solid-liquid 
case.  The  limitations  of  the  Bertoni  and  Tamir’s  formulation 
also  stem  from  the  fact  that  the  mathematical  analysis  be¬ 
comes  prohibitively  involved  as  the  system  to  be  considered 
becomes  more  complicated. 

Manuscript  received  January  28,  1980;  revised  March  28,  1980.  This 
work  was  sun  ’rted  by  the  Office  of  Naval  Research,  U.S.  Navy. 

The  autho  r  .ne  with  the  Physics  Department,  Georgetown  Univer¬ 
sity,  Washir  r  DC  20057. 
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Fig.  1.  Coordinate  system  and  typical  nonspecular  profiles  of  reflected 
and  transmitted  beams. 

Recently,  another  method  which  made  use  of  numerical 
integration  was  developed  by  Ngoc  [6]  to  describe  nonspec¬ 
ular  reflection  phenomena  for  a  liquid-solid  system.  This 
method  was  used  to  confirm  the  nonspecular  reflection  fea¬ 
tures  for  incidence  near  the  Rayleigh  critical  angle  and  to 
establish  nonspecular  reflectivity  for  incidence  near  the  longi¬ 
tudinal  critical  angle. 

The  present  paper  is  intended  to  apply  in  general  the  nu¬ 
merical  integration  method  to  investigate  the  nonspecular 
reflection  and  transmission  effects  for  a  layered  system.  Re¬ 
sults  for  incidence  near  critical  angles  of  a  liquid-solid  system 
will  be  discussed,  and  nonspecular  characteristics  of  a  liquid- 
solid-liquid  system  will  be  investigated  for  both  reflected  and 
transmitted  beams. 

II.  Theoretical  Framework 
A.  Spectral  Representation  of  a  Bounded  Beam 

Fig.  1  shows  an  incident  ultrasonic  beam  bounded  in  the 
(x,  z)  plane  and  uniform  in  thej  direction.  The  beam  of  an¬ 
gular  frequency  w  =  2rr/is  seen  to  be  incident  at  an  angle  0,. 
Its  intensity  profile  is  contained  by  an  effective  width  2w, 
which  is  projected  onto  the  first  interface  as  2w0 .  The  same 
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figure  also  illustrates  the  coordinate  system  for  a  liquid-solid  - 
liquid  structure.  The  same  coordinate  system  is  still  useful  for 
a  liquid-solid  system  when  the  second  interface  is  moved  to 
infinity. 

In  the  context  of  spectral  analysis,  a  bounded  beam  can  be 
represented  as  a  superposition  of  an  infinite  number  of  plane 
waves  which  have  different  amplitude  and  wave  vector.  It  is 
shown  [7]  by  using  Fourier  analysis  that  the  field  of  an  inci¬ 
dent  beam  can  be  uniquely  determined  at  any  point  ( x ,  z)  if 
its  field  distribution  is  known  in  any  plane.  Thus  a  spectral 
representation  of  a  well-defined  incident  beam  can  be  ex¬ 
pressed  by  the  following  Fourier  integral  transform  pair 


/*  kj  +  7T/  Wq 

Uinc(x,  z)  *  (2*)-1  I  exp  [iC**,  +  zkz)]  dkx 

■ikj-ir/w  o 


/wo 

t/j„c(x,0)  exp  (-  ixkx)  dx. 
w0 


(1) 

(2) 


In  (1 )  and  (2),  Uinc  (x,  0)  is  the  field  distribution  of  the  in¬ 
cident  beam  at  the  first  interface  where  z  =  0\kx  is  the x  com¬ 
ponent  of  the  wave  vector  k,  kz  is  the  z  component  of  k  de¬ 
fined  by  kz  =  ( k 2  -  kj)112 ;  the  wave  vector  kj  =  (co/d)  sin  0,- 
indicates  the  central  direction  of  propagation  of  the  incident 
waves;  V(kx)  is  interpreted  as  the  amplitude  of  the  constituent 
plane  waves.  The  time  dependence,  exp  (-  icof),  is  suppressed. 
In  this  paper,  the  incident  beam  is  taken  to  have  a  Gaussian 
distribution,  and  its  sound  field  at  the  first  interface  is  char¬ 
acterized  by 

^inc  (*,  0)  =  exp  [-  (x/w0 )2  +  ixkx  ]  (3) 


and  hence  from  (2) 

V(kx)  =  ff1/2w o  exp  [- (kx  -  kj)2  (w0/2)2 ] .  (4) 

The  physical  interpretation  of  (1)  is  that  the  incident  beam 
is  composed  of  an  infinite  number  of  plane  waves  having  the 
same  wavelength  but  incident  at  different  angles;  in  other 
words,  the  wave  vector  kx  is  perturbed  in  a  narrow  range 
about  kj.  The  principle  of  spectral  representation  can  be  ex¬ 
tended  to  describe  a  bounded  beam  reflected  from  or  trans¬ 
mitted  through  layered  media.  If  P(kx)  denotes  the  plane 
wave  reflection  or  transmission  coefficients  for  a  particular 
layered  structure,  then  the  associated  sound  field  can  be 
represented  by 

r  kj+*l  h»o 

U(x,z)  =  (2rr)-*  P(kx)  V(kx) 

Jkj-n/w o 

■  exp  [/ (jcAcjj  +  zkz))  dkx.  (5) 


B.  Plane  Wave  Reflection  or  Transmission  Coefficients 
It  is  evident  from  (5)  that  the  relevant  plane  wave  coefficient 
P(kx)  strongly  influences  the  beam  profile.  Indeed,  (5)  indi¬ 
cates  that  the  resulting  profile  is  constructed  by  the  interfer¬ 
ence  of  the  individual  plane  waves,  which  have  amplitude 
P(kx)  V(kx)  and  undergo  a  phase  shift  upon  reflection  or 


transmission.  In  this  paper,  one  is  interested  in  the  profile  of 
the  sound  beam  reflected  from  a  liquid-solid  interface  and 
those  of  both  reflected  and  transmitted  beams  observed  in  the 
liquid-solid-liquid  system. 

The  amplitude  plane  wave  coefficients  P(kx )  for  the  cases  of 
interest  are  given  [8]  by 

RLs(kx)  =  N/D  (6) 

the  reflection  coefficient  for  a  liquid-solid  infinite  half  space 
system;  and 

^lsl  (kx)=NK/FaFz  (2) 

and 

TtSL(kx)=NTIFaFs  (8) 

for  reflection  and  transmission  for  a  liquid-solid-liquid  system. 
The  symbols  used  in  (6)-(8)  are  defined  as  follows: 

N  =  (kj  -  2 k2x)2  +  4 k2xKsKd  -  pkjKJK 

D  =  (kj  -  2kj)2  +  4 klKsKd  +  pkjKJK 

Nr  =(kj  -  2ki)4  +  16 k'KjK2,  -  p2kssK2dIK 2 

+  8 (kj  -  2k2x)2k2xKsKd(\  -  cos  Pcos  Q)/(sin  Psin  Q) 

Nt  =  (lipkj  KJK)  [(*|  -  2fc3)J/sin  P  +  4k2xKsKdlsm  Q ] 

Fa  =  (kj  -  2 kj)2  (1  -  cos  P)/sw  P 

+  4k2xKsKd(l  -  cos(2)/sinQ  +  ipkjKdIK 

Fs  -  (kj  -  2 kj)2  (1  +  cos  P)l sin  P 

+  4k2xKsKd(\  +  cos  0)/sin  Q  -  ipkj  Kd/K 

where 

—  '  I  " 

P-PIP 

k  =  (w/u)  (1  +  ial2it) 
kd  =(w/ud)(l  +iad/2n) 

*s  =  (w/ u,)(l  +iaj  2ir) 

k  =  (**  -  k\y >2 

K*  k2d-k2xyl 2 

Ks=(kj-k2xyi2 
P  =  dKd 
Q=dKs. 

The  quantities  a,ad,  and  as  denote  the  loss  parameters  defined 
as  attenuation  per  wavelength  for  the  sound  waves  in  liquid 
and  for  the  longitudinal  and  shear  waves  in  solid, respectively. 

It  can  be  readily  shown  that  R\,s(kx)  is  a  function  of  inci¬ 
dent  angle  and  frequency,  while  R LSl  (** )  and  TLSL  (kx)  are 
functions  of  incident  angle  and  the  product  fd.  The  frequency 
dependence  is  present  because  the  loss  parameters  depend  on 
the  sound  frequency. 

In  general,  the  coefficients  P(kx)  are  complex  functions 
possessing  pole-zero  pairs  which  completely  determine  their 
analytical  behavior  in  the  complex  plane  [8] .  These  poles  and 
zeros  can  be  arranged  into  pairs  because  each  consists  of  one 
pole  and  one  zero,  both  having  the  same  real  part  in  the  loss¬ 
less  case.  For  a  liquid-solid-liquid  system,  the  common  real 
part  of  a  pole-zero  pair  is  associated  with  an  angle  of  incidence 
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6,0 

Fig.  2.  Moduli  for  lossless  water-brass-water  system  with  fd  -  2.4  km/s. 
(a)  rLSl(*jc)-  (b)  TL$i(kx). 


20  40  60 

0,0 

Fig.  3.  Phases  for  lossless  water-brass-water  system  with  fd  =  2.4  km/s. 
WRLSL<fcxl-  (b)rLSL^x>- 

where  a  plate  mode  can  be  excited.  This  association  is  ex¬ 
pressed  by  the  following  relationship 

l*pl/l*l  =  l*ol/l*l  =  sin0p  (9) 

where  \k  |  =  w/o,  0p  is  the  plate-mode  angle  of  incidence,  and 
kp  and  k0  are  pole  and  zero  of  the  same  pair.  At  a  plate-mode 
angle  of  incidence,  |/?LSL  I  =  0  and  I  T^sl  I  =  1  when  the 
media  are  considered  lossless,  and  the  phase  curves  of  /?LSL 
and  Tlsl  exhibit  a  point  of  inflection.  Fig.  2  shows  |f?LSL  I 
and  I TLSL  |  as  functions  of  incident  angle  for  the  lossless 
water-brass-water  system  with  fd  =  2.4  km/s.  The  relative 
phase  of  Klsl  and  Tlsl  f°r  the  same  system  is  plotted  versus 
angle  of  incidence  in  Fig.  3. 


20  40  60  80  100 

2fw(km/s  1 

Fig.  4.  Half-width  Wv  of  Pfk*)  plotted  against  product  frequency  times 
beamwidth.  Value  used  for  9,  is  30°. 


C.  Functional  Characteristics  of  V(kx) 

Equation  (5)  also  indicates  that  the  resulting  profiles  are  de¬ 
termined  by  the  functional  behavior  of  V(/c*),  which  is  the 
Fourier  transform  of  the  incident  intensity  distribution.  It  is 
seen  from  (4)  that  )  is  a  Gaussian  function  whose  shape 
depends  on  beamwidth,  sound  frequency,  and  angle  of  inci¬ 
dence.  With  the  values  of  beamwidth  being  a  few  centimeters 
and  those  of  frequency  ranging  from  1  MHz  to  10  M  •  -,  the 
magnitude  of  V(kx)  is  significant  only  in  a  narrow  >a  ige  about 
kj  and  almost  independent  of  the  values  of  inciden.  angle.  It 
is  reasonable  to  represent  the  significant  range  of  the  Gaussian 
function  F(kx)  by  its  half-width.  This  half-width  is  denoted 
by  Wv  and  can  be  written  in  units  of  degrees  as 

W„  =  sin"1  (sin  0,  +  0.53  W*)  -  sin'1  (sin  0,  -  0.53  Wk) 

(10) 

where 

Wk  =  v  cos  0,/2/w.  (11) 

Fig.  4  illustrates  the  functional  dependence  of  Wv  versus  the 
product  frequency  times  beamwidth  for  0,-  =  30°.  Considering 
the  known  forms  of  P{kx)  and  V(kx),  one  can  conclude  that 
the  calculated  profile  is  a  function  of  0(,/,  and  2 fw  in  the 
liquid-solid  case,  and  of  0,-,/d,  2 fw  for  a  liquid-solid-liquid 
system. 

D.  Alternative  Methods  of  Calculating  Beam  Profiles 

To  simplify  the  mathematical  involvement  and  to  avoid  the 
complication  that  a  beam  profile  may  change  as  the  beam 
travels  away  from  the  interfaces  [3] ,  the  beam  profiles  for 
the  cases  of  interest  are  all  calculated  at  the  interfaces,  that  is, 
at  either  z  =  0  or  z  =  -d.  Under  such  circumstances,  the  sound 
field  expressed  by  (5)  becomes 

rkj  +  */w0 

l/LS(x,  0)  =  (2n)-1  Rts(kx)V(kx) 

Jkj-n /w0 

•  exp(ikxx)dkx  (12) 

rkj*  ir/w0 

U?si  (x.  0)  =  (2rr)" 1  RLSL  (kx)  V(kx) 

Jki-*lw o 


■  exp  ( ikxx)dkx 


(13) 
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rkj-ml  tv0 

U[SL(x,-d)  =  (2n)-'  I  TLSL(kx)  V(kx) 

JKj  -  <r/w0 

•  exp  kzd)]  dkx<  (14) 

where  t/LS(x,  0)  is  the  reflected  field  for  a  liquid-solid  system, 
and  (JlSL(x ,  0)  and  U\^sh(x,- d)  are  the  reflected  and  trans¬ 
mitted  fields  for  a  liquid-solid-liquid  system. 

Two  methods  are  known  for  the  evaluation  of  (12)-(14). 
Bertoni  and  Tamir  [1  ] ,  Pitts  [8] ,  and  Ng  [5]  chose  to  evalu¬ 
ate  these  integrals  analytically  after  replacing  the  exact  but 
complicated  espressions  of  P(kx )  with  approximate  ones  con¬ 
structed  from  the  knowledge  of  the  location  of  one  or  several 
pole-zero  pairs  they  considered  significant.  This  method  will 
be  referred  to  as  the  Bertoni  and  Tamir  method.  The  other  ap¬ 
proach  used  by  Ngoc  [6|  is  to  evaluate  integrals  numerically. 
The  numerical  integration  method  has  been  shown  to  be  suc¬ 
cessful  for  some  simple  cases  [9] ,  [10] . 

The  Bertoni  and  Tamir  method  essentially  proceeds  as 
follows:  1)  determine  the  locations  of  pole-zero  pairs  of  P(kx)\ 

2)  formulate  a  simple  approximate  expression  for  P(kx);  and 

3)  carry  out  the  still  rather  involved  integration  analytically. 
The  numerical  integration  method  proceeds  in  a  more  direct 
manner,  numerically  evaluating  the  integrals  on  the  basis  of 
the  Simpson’s  integration  algorithm  thus  omitting  the  mode- 
identifying  step. 

The  numerical  integration  method  is  found  to  be  efficient 
and  adaptable  to  some  practical  considerations  such  as  beam 
spreading  effect  and  non-Gaussian  incident  profiles.  In  addi¬ 
tion,  the  accuracy  of  the  numerical  integration  method  is  not 
restricted  to  angles  of  incidence  corresponding  to  resonant 
modes  as  is  the  case  of  the  Bertoni  and  Tamir  method.  All 
calculated  results  presented  below  are  obtained  by  the  nu¬ 
merical  integration  method. 

III.  Calculated  Proeiles  eor  a  Liquid-Solid 
System 

It  has  been  established  that  for  a  liquid-solid  system  [6] 
nonspecular  phenomena  are  prominent  at  such  an  angle  of  in¬ 
cidence  that  either  or  both  of  the  modulus  and  phase  of 
R^s(kx)  vary  strongly.  Such  behavior  is  found  near  the  longi¬ 
tudinal  critical  angle  0d  for  a  water-Plexiglas  interface  and 
near  the  Rayleigh  critical  angle  0r  for  various  liquid-solid  in¬ 
terfaces.  The  reflected  beam  profile  is  a  function  of  the  ab¬ 
sorption  as  well  as  of  0hf,  and  2 fw.  In  the  following  calcula¬ 
tions  w  and  0j  are  varied, but  / remains  fixed.  Therefore, 
profile  changes  as  w  is  varied  can  be  associated  with  V(kx), 
while  changes  resulting  from  variation  of  0,  can  be  attributed 
to  R(kx). 

A.  Incidence  near  Longitudinal  Critical  Angle 

Fig.  5  shows  the  variation  of  the  modulus  and  phase  of 
Ris(kx)  near  6d  for  a  water-Plexiglas  system  at  different  fre¬ 
quency  values.  Both  sets  of  curves  exhibit  strong  variations 
in  the  angular  range  under  consideration.  Profiles  of  the  re¬ 
flected  beam  are  calculated  by  the  numerical  integration 


Fig.  S.  Absorptive  near  longitudinal  critical  angle  for  water- 

Plexiglas  interface,  (a)  Modulus,  (b)  Phase. 


X/w„ 

Fig.  6.  Water-Plexiglas  if  -  2  MHz):  calculated  reflected  beam  profiles, 
(a)  8j  =  Bd  as  beamwidth  is  varied  (4:  1 2.7  mm.  H  19,05  mm,  and  C: 
25.4  mm),  (b)  Beamwidth  is  19.05  mm  as  0,  increases  beyond  0d 
(4:  Sd  +  1°,B:  8d  +  2 ,C  0d+  1,  and  D:  6d). 

method  for/=  2  MHz  when  tv  and  0,-  are  changed.  Non¬ 
specular  reflectivity  is  manifest  only  in  the  form  of  a  lateral 
beam  displacement,  the  magnitude  of  which  varies  signifi¬ 
cantly  with  both  w  and  0,.  These  results  are  presented  in 
Fig.  6.  It  is  noted  that  the  maximum  beam  displacement  does 
not  occur  exactly  at  6d  but  approximately  at  6d  +  1°. 

B.  Incidence  near  Rayleigh  Critical  Angle 
In  Fig.  7,  the  modulus  of  R(kx)  is  plotted  versus  0,  near  8r 
for  a  water-stainless  steel  system.  The  phase  curves  are  not 
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. . .  lossless 
_  absorptive 


Fig.  7.  Modulus  of  absorptive  R ls(*x)  near  the  Rayleigh  critical  angle 
for  water-stainless  steel  interface. 


X/w» 

Fig.  8.  Water-stainless  steel  (/=  2  MHz):  calculated  reflected  beam 
profiles,  (a)  8,-  =  6r  as  beamwidth  is  varied  04 :  12mm,  B-  20  mm, 
and  C  28  mm),  (b)  Beamwidth  is  20  mm  as  Br  is  varied  { A :  Br, 

B  Br  +  0.5°,  and  C  6r  +  1°). 


X/w. 

Fig.  9.  Water-brass-water  (fd  =  2.4  km/s)  beam  profiles,  (a)  Reflected, 
(b)  Transmitted.  Calculated  for  =  2.4°  at  various  plate-mode 
angles  of  incidence  (.4:  5 °,B:  24°,  C  45°,  and  D.  50°). 


IV.  Calculated  Profiles  for  a 
Liquid-Solid-Liquid  System 
The  numerical  integration  method  is  now  applied  to  study  a 
more  complicated  system  of  a  solid  plate  immersed  in  a  liquid. 
The  integrals  to  be  evaluated  are  (13)  and  (14),  which  describe 
the  intensity  profiles  of  the  reflected  and  transmitted  beams, 
respectively.  As  mentioned  in  Section  II,  the  physical  parame¬ 
ters  of  interest  in  this  case  are  8hfd,  and  2 fw.  To  separate  the 
influence  of  V{kx )  and  P(kx ),  f  is  kept  fixed  at  2  MHz,  and  w 
and  9j  are  varied. 

It  is  interesting  to  investigate  the  profiles  for  different  ranges 
of  incident  angles:  those  corresponding  to  plate  modes  and 
those  falling  between  two  adjacent  modes.  The  former  will  be 
referred  to  as  plate-mode  incidence  and  the  latter  as  between- 
mode  incidence.  For  this  reason,  calculations  for  a  water- 
brass-water  system  are  done  throughout  for  two  values  of  fd, 
2.4,  and  7.0  km/s,  since  the  plane  wave  coefficients  exhibit 
different  characteristics  in  terms  of  number  of  modes  and 
spacings  between  modes  for  these  values  of  fd.  In  the  follow¬ 
ing,  between-mode  incidence  is  used  to  indicate  the  angle  of 
incidence  exactly  between  two  adjacent  plate  modes,  whose 
angular  separation  is  denoted  by  S. 


shown  since  they  do  not  deviate  significantly  from  the  lossless 
case.  Again  the  reflected  beam  profile  is  calculated  for  differ¬ 
ent  values  of  w  and  8h  holding/fixed  at  2  MHz.  The  non- 
specular  characteristics  consist  of  an  intensity  minimum,  a 
beam  displacement,  and  a  trailing  sound  field.  As  the  beam- 
width  becomes  larger,  the  prominent  change  in  the  profile  is 
a  shift  of  the  whole  beam  to  the  right.  For  variation  of  9h 
the  most  noticeable  profile  change  is  seen  at  the  intensity  min¬ 
imum  which  gradually  gains  in  intensity  and  finally  disappears 
approximately  at  8r  +  1°.  These  results  are  presented  in  Fig.  8. 


A.  Plate-Mode  Incidence 

The  profiles  are  first  calculated  at  the  plate-mode  angles  of 
incidence  using/rf  =  2.4  km/s.  The  value  of  beamwidth  is 
taken  to  be  19.05  mm,  which  gives  2 fw  =  38.1  km/s.  It  is 
estimated  from  Fig.  4  that  Wv  is  equal  to  2.4°.  Thus  the  beam 
profiles  calculated  with  these  parametric  values  would  show 
only  the  functional  behavior  of  P(kx )  near  a  plate  mode,  be¬ 
cause  S  is  substantially  larger  than  IV,,  (see  Figs.  2  and  3).  The 
reflected  profiles  plotted  in  Fig.  9  show  features  similar  to 
those  present  in  the  beam  reflected  near  8r  from  a  liquid-solic 
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e,o 

Fig.  10.  Moduli  for  lossless  water-brass-water  system  with/d  =  7.0 
km/s.  (a)  TLSL(kx).  (b)  /tLSL<**)- 

interface  (compare  Fig.  8).  The  same  figure  also  shows  the 
profiles  calculated  for  the  transmitted  beam,  displaying  only 
a  simple  beam  shift  for  all  modes.  The  absence  of  an  intensity 
minimum  in  the  transmitted  profiles  can  probably  be  attrib¬ 
uted  to  the  fact  that  the  phase  of  T l.sl(^x  )  does  not  exhibit 
a  360°  change  as  ^lsl(^jc)  normally  does  at  plate-mode  in¬ 
cidence.  It  is  also  noted  that  the  intensity  of  the  transmitted 
profiles  increases  with  the  width  of  the  plate-mode  peaks  of 
the  modulus  curves. 

B.  Between-Mode  Incidence 

In  this  angular  range  fd  =  7.0  km/s  is  used  so  that  S  varies 
noticeably  for  the  same  fd  (Fig.  10).  The  chosen  value  of  Wv 
is  still  2  .4°.  Both  reflected  and  transmitted  profiles  are  cal¬ 
culated  at  three  between-mode  angles  of  incidence  selected  to 
be  15.6°,  24.7°,  and  41 .4°,  whose  values  of  S  are  0.75°,  5.35°, 
and  3.6°,  respectively.  Fig.  1 1  shows  the  resultant  profiles.  It 
is  observed  that  the  profiles  A  display  features  characteristic 
of  a  plate-mode  incidence  profile,  those  labeled  B  have  no 
noticeable  nonspecular  features,  and  the  profiles  C  show  an 
“interference”  effect  caused  by  the  two  adjacent  modes. 
Comparison  of  the  values  of  Wu  and  S  leads  to  the  following 
proposition :  interference  effect  is  most  probable  when  Wu  < 

S  <  2WV\ nonspecular  characteristics  are  not  likely  when 
S  >  2  ;  and  when  S  «  Wv  a  nonspecular  profile  produced 

in  the  neighborhood  of  two  adjacent  modes  is  similar  to  that 
resulting  for  plate-mode  incidence. 

C.  Variation  of  Incident  Angle 

The  basic  characteristics  of  the  reflected  and  transmitted 
profiles  have  been  established  for  plate-mode  and  between- 


Fig.  11.  Water-brass-water  ( fd  ~  7.0  km/s)  beam  profiles,  (a)  Re¬ 
flected.  (b)  Transmitted.  Calculated  for  Wv  =  2.4°  at  various  be¬ 
tween-mode  angles  of  incidence  (A'-  15.6 °,  B:  24.7°,  and  0  41.4°). 


X/  Wo 

Fig.  12.  Water-brass-water  (fd  -  2.4  km/s)  beam  profiles,  (a)  Re¬ 
flected.  (b)  Transmitted.  Calculated  for  Wv  =  2.4°  and  9,  being 
varied  about  plate-mode  incident  angle  9^  -  24°  (A :  24°,  B.  25°.  and 
C:  26°). 

mode  types  of  incidence.  In  this  subsection,  these  are  investi¬ 
gated  to  see  how  they  change  with  the  angle  of  incidence. 
Among  those  of  Figs.  9  and  1 1 ,  the  profile  for  0,  =  24°  is 
chosen  for  the  case  of  plate -mode  incidence,  and  the  one 
which  exhibits  the  interference  effect  is  chosen  for  investiga¬ 
tion  of  the  between-mode  case.  It  is  evident  from  Figs.  1 2  and 
13  that,  when  0,  is  varied  from  its  initial  value,  the  nonspec¬ 
ular  features  become  less  pronounced  and  the  nonspecular 
profiles  are  gradually  transformed  into  specular  ones.  Fig.  13 
in  particular  indicates  that  the  interference  peaks  are  less 
prominent  as  the  incident  angle  deviates  from  the  between- 
mode  incidence  of  0<  =  41 .4°. 
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Fig.  13.  Water-brass-water  (fd  =  7.0  km/s)  beam  profiles,  (a)  Re¬ 
flected.  (b)  Transmitted.  Calculated  for  Wv  =  2.4°  and  0;  being 
varied  about  between-mode  incident  angle  41.4°  (A:  41.4 41.1°, 
and  C  40.8°). 


X/w. 

Fig.  15.  Water-brass- water  (fd  -  7.0  km/s)  beam  profiles,  (a)  Re¬ 
flected.  (b)  Transmitted.  Calculated  at  between-mode  incidence 
6 j  =  41.4°  for  various  values  of  Wv  (A  :  2.4°.  B:  3.55°,  and 
C:  7.2°). 


X/w. 

Fig.  14.  Water-brass-water  (fd  =  2.4  km/s)  beam  profiles,  (a)  Re¬ 
flected.  (b)  Transmitted.  Calculated  at  plate-mode  incidence  = 
45°  for  various  values  of  (A :  2.4°,  B.  3.55°,  and  C:  7.2°). 

D.  Variation  of  Beamwidth 

Finally,  the  profile  changes  are  examined  when  the  beam- 
width  is  varied.  Variation  of  the  beamwidth  amounts  to  a 
change  in  the  half-width  of  V(kx).  The  beamwidth  values 
used  in  the  following  calculations  are  19.05  mm,  12.7  mm, 
and  6.35  mm.  Therefore,  the  values  of  W„  derived  from  Fig 
4  are  2.4°,  3.55°,  and  7.2°,  respectively.  Again  calculations 


are  done  for  both  types  of  incidence,  0,  =  45°  chosen  for  the 
plate-mode  case  and  0,-  =  41 .4°  for  the  between-mode  case.  In 
the  plate-mode  case  (Fig.  14),  the  profiles  are  qualitatively  not 
changed  when  is  equal  to  2.4°  and  3.55°.  For  Wv  =  7.2°, 
both  reflected  and  transmitted  profiles  have  an  additional  peak 
caused  by  the  influence  of  the  nearest  mode  (0,  =  50°,  i.e.,  5° 
away),  which  is  partially  overlapping  the  mode  under  consider¬ 
ation.  The  results  for  between-mode  incidence  are  illustrated 
in  Fig.  15.  The  interference  effect  is  still  present  for  the  low 
values  of  Wv.  When  =  7.2°,  both  reflected  and  transmitted 
profiles  show  some  plate-mode  incidence  character  since  the 
half-width  Wu  is  now  sufficiently  larger  than  S  which  is  equal 
to  3.6°  in  this  case. 


V.  Conclusion 

The  numerical  method  described  above,  when  applied  to  cal¬ 
culation  of  reflection  from  a  liquid-solid  interface,  shows  that 
nonspecular  reflection  of  a  Gaussian  incident  beam  can  occur 
at  both  the  longitudinal  and  Rayleigh  critical  angles.  The 
strength  and  observability  of  the  effect  depends  on  the  fre¬ 
quency,  the  beamwidth,  and  the  mechanical  properties  of  the 
liquid-solid  combination  chosen. 

The  same  method  applied  to  a  liquid-solid-liquid  system 
yields  information  about  nonspecular  reflectivity  and  trans- 
mittivity.  It  is  concluded  from  the  results  presented  that  0,- 
and  fd  are  not  the  only  parameters  of  importance  but  that  a 
new  parameter,  2fw,  the  product  frequency  times  beamwidth. 
determines  whether  nonspecular  reflection  and  transmission 
can  be  expected  to  occur  at  incident  angles  other  than  plate¬ 
mode  angles  of  incidence. 
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Detection  of  local  inhomogeneities  in  solids 
by  Rayleigh  angle  reflection 

W.G.  MAYER 

A  schlieren  method  is  employed  to  observe  changes  in  the  reflection  profile  of  an 
ultrasonic  beam  incident  at  the  Rayleigh  angle.  These  changes  are  used  to  determine 
the  location  of  small  elastic  inhomogeneities  of  a  flat  solid  sample. 


Introduction 

Ultrasonic  methods  have  been  used  for  many  years  to 
determine  not  only  the  elastic  properties  of  solids  but  also 
changes  of  these  properties  and  inhomogeneities  which 
influence  propagation  velocities  in  the  solid.  As  early  as 
1946,  Firestone  and  Frederick1  used  ultrasonic  velocity 
measurements  to  find  anisotropies  in  rolled  aluminium. 
Changes  in  ultrasonic  velocities  can  be  caused  by  a 
variety  of  anisotropies,  among  them  changes  in  Young’s 
modulus5  which  may  be  measurable  as  variations  in  the 
shear  wave  velocity.3  Longitudinal  wave  velocity  changes 
can  be  caused  by  preferred  orientations,  as  for  instance 
in  a  drawn  zinc  bar  where  the  velocity  was  found  to  vary4 
between  3540  m  s'1  and  4400  m  s'1 . 

In  addition  to  bulk  wave  velocity  measurements  one  can 
use  surface  waves  to  determine  changes  in  bulk  properties5 
or  texture.6  Lamb  waves  can  also  be  used1  for  the 
determination  of  irregularities  in  thin  extended  solids. 

Generally,  measurements  made  using  these  techniques 
reveal  changes  in  sound  velocities  caused  by  anisotropies 
that  exist  either  throughout  the  entire  length  of  the  sound 
path  or  within  a  portion  of  the  measurement  distance. 
However,  a  small  isolated  irregularity  may  contribute  very 
little  to  the  overall  change  in  velocity  determined  over  a 
relatively  long  distance  which  includes  the  irregularity. 

The  method  described  below  samples  elastic  properties 
of  the  solid  and  reveals  small  local  changes  in  the  ultrasonic 
velocity  which,  if  a  long  path  were  used,  could  easily  escape 
delection. 

Theoretical  background 

The  method  described  is  based  on  the  fact  that  a  bounded 
ultrasonic  beam  incident  from  a  liquid  on  to  a  flat  solid 
may  be  non-specularly  reflected  if  the  angle  of  incidence 
is  equal  to  the  Rayleigh  angle.  One  of  these  non-specular 
reflection  phenomena,  a  lateral  beam  displacement,  was 
first  observed  by  Schoch."  Whether  the  reflected  beam  is 
laterally  displaced  or  whether  an  intensity  profile  change 
occurs,  that  is.  the  reflected  beam  appears  to  split  into  two 
distinct  portions  separated  by  an  intensity  null  region,9 
depends  strongly  on  the  beam  width,  the  frequency,  and 

The  author  is  prpsently  at  the  Labor  Kur77pi!physik.  Universitat 
des  Saariandes,  6600  Saarbrucken,  West  Germany.  H  $  permanent 
address  ts  the  Physics  Department,  Gecrcetown  University, 
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on  how  closely  the  incidence  angle  approaches  the 
Rayleigh  angle.10 

These  changes  in  the  reflected  beam  profile  can  be  calculated 
if  one  considers  the  reflection  coefficient  for  a  bounded  beam 
rather  than  a  single  plane  wave.  Representing11  the  incident 
beam  as  a  summation  of  an  infinite  number  of  plane  waves, 
all  of  the  same  frequency  but  incident  at  the  boundary  at 
different  angles,  one  can  construct  the  reflected  sound  field 
at  the  boundary  for  any  point  along  the  interface,  that  is, 
the  x-direction.  The  amplitude  distribution  of  the  reflected 
beam  is  given12  by 


Unix)  =  (1/:tt)  j  R(kx)  V(kx)  exp(ix^)  dkx  (1) 


where  R(kx)  is  the  plane  wave  amplitude  reflection 
coefficient.  kx  is  the  projection  of  the  incident  wave  vector 
on  the  x-axis.  and  l’(kA  ).  the  Fourier  transform  of  the 
incident  particle  displacement,  depends  on  the  incident 
field  profile.  U^fx), 


•  (**)-  I  tWIxjexpL-  irkx)  dx.  (21 


Introducing  attenuation,  (he  various  terms  in  the  above 
expressions  can  be  transformed  in  such  a  fashion  that  the 
reflected  beam  profile  can  be  calculated  bs  numerical 
integration.13  Numerical  evaluations  for  different  flat 
solids  immersed  in  water  show  that  a  Gaussian  beam 
incident  at  the  Rayleigh  angle  is  reflected  with  an  intensity 
profile  that  is  quite  different  from  a  Gaussian  distribution. 
The  resulting  profile  depends  on  beam  width  and 
frequency  and.  most  important  for  the  present  purpose,  on 
whether  the  incidence  is  exactly  at  the  Rayleigh  angle  or 
at  a  slightly  different  angle.  This  latter  dependence  is 
illustrated  for  the  case  of  a  water-stainless  steel  boundary 
as  shown  in  Fig.  1 ,  which  is  adapted  from  Ngoc.10  The 
centre  of  the  incident  Gaussian  beam,  not  shown  in  Fig.  1 . 
is  incident  at  point  0  on  the  x-axis.  The  graph  indicates 
that  the  reflected  beam  profile  is  no  longer  Gaussian  and 
chances  markedly  when  the  angle  of  incidence  is  changed 
slightly  from  6R.  the  Rayleigh  angle. 
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Fig.  1  Calculated  profile  of  reflected  beam  at  a  water  stainless 
steel  boundary  for  incident  angles  of:  A,  0R;  B,  +  0,5“ ; 

C,  +  1c.  Beam  width  is  20  mm,  frequency  2  MHz,  and  incident 
beam  (not  shown)  is  centred  at  X  !W0  -  0,  where  WQ  is  a  measure 
of  beam  width.  (Aftet  Ngoc10) 

Experimental  method  and  results 

The  curves  in  Fig.  1  arc  given  as  a  function  of  Oj,  the  angle 
of  incidence,  for  0j  =  0R  to  0 ,  =  +  1  °.  The  Rayleigh  angle 

for  stainless  steel  is  30.6°.  determined  by  the  Rayleigh 
wave  velocity.  vR.  and  the  sound  velocity  in  water,  v. 
according  to 

sin  OJr  =  l/rR.  (3) 

With  v  =  1490  m  s'1  and  0,  =  30.6°  one  obtains  rR  = 

2930  ms*1. 

If  one  immerses  a  flat  sample  of  stainless  steel  in  water,  a 
2  MHz.  20  mm  wide.  Gaussian  ultrasonic  beam  impinging 
at  0,  =  Or  will  be  reflected  according  to  curve  A  in  Fig.  1 . 
Such  a  profile  can  be  observed  visually  with  a  standard 
Schlieren  arrangement.  If  one  now  moves  the  sample 
parallel  to  its  surface  the  profile  will  not  change  if  i’R 
remains  the  same  at  all  locations  subjected  to  the  incident 
beam.  If.  however,  the  beam  profile  observed  at  some 
location  deviates  from  the  previously  observed  shape, 
the  value  of  rR  in  (3)  must  be  different  at  that  particular 
location  since  neither  0j  nor  r  have  changed. 

One  can  easily  observe  on  a  homogeneous  sample  changes 
in  the  reflection  profile  caused  by  a  change  in  the  angle 
of  incidence.  A0j,  of  0.25°  (see  Fig.  1 )  if  one  were  to  change 
the  angle  of  incidence  without  moving  the  sample.  Con¬ 
versely,  one  can  also  notice  quite  readily  changes  in  the 
profile  caused  by  a  change  in  the  Rayleigh  wave  velocity. 
ArR.  if  one  does  not  change  0,  but  moves  the  sample, 
provided  the  sample  is  locally  inhomogeneous.  An 
example  of  this  situation  is  shown  in  Fig.  2.  Here  a  10  mm 
wide.  2  MHz  beam  impinges  on  a  flat  (to  at  least  0.05°) 
thick  brass  sample.  The  angle  of  incidence  is  the  Rayleigh 
angle  and  the  resulting  reflection  profile  is  shown  in 
Fig,  2a.  Moving  the  sample  parallel  to  its  surface  docs  not 
change  the  appearance  of  the  reflected  beam  profile  until 
one  reaches  an  area  on  the  sample  where  the  elastic 
properties  of  the  sample  change,  giving  rise  to  a  different 
non  specular  reflection  profile,  shown  in  Fig.  2b.  The 
original  appearance  can  be  restored  by  a  slight  adjustment 
of  0,.  In  the  case  of  this  particular  sample  a  decrease  of 
0.75°  in  0 ,  was  required. 

Substituting  this  adjusted  angle.  0,  =  Or  0.75°,  into  (3) 
and  comparing  the  resulting  value  of  rR  with  that  one 
obtains  for  the  rest  of  the  sample  one  finds  that  the  local 


Fig.  2  Schlieren  photographs  of  2  MHz  beam  reflected  at  a 
brass-water  flat  interface  with:  a  -  Rayleigh  angle  incidence;  and 
b  -  same  angle  of  incidence  but  at  section  of  sample  with  slightly 
different  elastic  properties 

value  of  rR  is  about  24  m  s"1  higher  than  elsewhere  on  the 
sample.  Thus,  this  method  is  sensitive  to  small,  localized 
changes  in  the  Rayleigh  wave  velocity. 

Discussion 

The  method  described  above  can  be  used  to  find  local 
Rayleigh  wave  velocity  fluctuations  that  can  be  as  small 
as  a  fraction  of  one  percent  of  rR.  Assuming  that  a  profile 
change  is  detectable  that  corresponds  to  a  0.25°  deviation 
from  0H,  the  resulting  detectability  of  a  change  in  »>R 
depends  on  the  magnitude  of  the  unperturbed  Ray  leigh 
wave  velocity.  Fig.  3  shows  the  amount  of  ArR  for  A0j  of 
0.25°.  0.5°.  and  1°  as  a  function  of  the  overall  i’R  and  the 
corresponding  0R.  The  curves  shown  are  for  a  solid-water 
boundary,  calculated  from  (3)  for  positive  A0;.  The 
corresponding  ArK  for  negative  A0,  are  shown  as  points. 

Although  this  method  enables  one  to  detect  small,  local 
changes  in  the  elastic  properties  one  should  not  expect  to 
be  able  also  to  determine  the  cause  of  these  velocity 
anisotropies.  Considering  that  rR  depends  on  the  value 
of  the  shear  wave  velocity  and  on  Poisson’s  ratio,  which 
itself  depends  on  the  shcat  wave  and  the  longitudinal  wave 
velocities,  any  observed  ArR  may  be  caused  by  local 
fluctuations  in  these  quantities.  Moreover,  since  these  bulk 
wave  velocities  depend  on  the  elastic  moduli  and  the 
density  of  the  solid,  an  observed  ArR  may  also  be  caused 
by  anisotropies  of  these  latter  parameters. 

Thus  the  method  described  here  can  be  primarily  useful 
in  the  location  of  small  localized  inbomogeneities  of  the 
solid  while  the  determination  of  the  exact  causes  can  then 
he  attempted  wilh  the  appropriate  techniques  applied  at 
the  location  in  question. 
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Fig.  3  Changes  in  detectable  with  A,  0.25c;  6 ,  0.5°; 

C,  1c  as  a  function  of  v R  or 
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